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Abstract. The propagation of the heavy quarks produced in relativistic nucleus-nucleus collisions at RHIC 
and LHC is studied within the framework of Langevin dynamics in the background of an expanding 
deconfined medium described by ideal and viscous hydrodynamics. The transport coefficients entering into 
the relativistic Langevin equation are evaluated by matching the hard-thermal-loop result for soft collisions 
with a perturbative QCD calculation for hard scatterings. The heavy-quark spectra thus obtained are 
employed to compute the differential cross sections, the nuclear modification factors Raa and the elliptic 
flow coefficients V2 of electrons from heavy-flavour decay. 



1 Introduction 

The aim of the ongoing heavy ion programme at RHIC 
(with Au-Au collisions at v^nn '^P -^OO GeV) and of 
the recently started programme at LHC (with the first 
Pb-Pb collisions at -s/snn = 2.76 TeV) is the creation of 
a deconfined system of quarks and gluons, with a lifetime 
sufficiently long (a few fm) to allow measurable signals to 
show up. The challenge is not only to establish the forma- 
tion of the deconfined phase, but also to study the proper- 
ties of the produced medium. Among the most interesting 
observables investigated at RHIC, perhaps, there are the 
quenching of high-p^ hadron spectra, expressed through 
the nuclear modification factor RAAiPr) [I1I2]; and the az- 
imuthal anisotropy of particle emission in the transverse 
plane, quantified by the elliptic-flow coefficient V2 [3,4 . 
The asymptotic value Raa{pt) ~ 0.2 found at high pr 
indicates that the matter created in the heavy ion colli- 
sions is very opaque, leading to a sizable energy loss of 
the few high-momentum partons produced in the initial 
instants of the collisions by hard processes. The latter 
are describable in the framework of perturbative Quan- 
tum Chromodynamics (pQCD). Furthermore, the elliptic 
flow observed in semi-central collisions seems to find a nat- 
ural explanation within a hydrodynamical description of 
the medium, displaying an expansion driven by the pres- 
sure gradient: the latter, resulting larger along the reaction 
plane, leads to the azimuthal anisotropy of the transverse- 
momentum spectra. Notably, the large value of V2 found 
in the data, besides supporting the hydrodynamical be- 
havior of the created matter (entailing a mean free path 
Amfp much smaller then the system size L), requires also 
a quite short thermalization time (to<1 fm). Raa and V2 



have recently been measured also by the ALICE detector 
at LHC for unidentified particles [SUB]. 

These two observables (i?AA an V2) have been initially 
studied for the case of light-flavour hadrons (pions, kaons, 
protons, ...), arising from the hadronization of gluons and 
light quarks. In this context, the quenching of hadron spec- 
tra has been attributed to the energy-loss of hard partons 
due to the radiation of soft gluons (uj ^ £'hard) occurring 
in the presence of scattering centers (the plasma particles) . 
Different theoretical approaches have been proposed [71[S1 
I^ITU] , all stressing the importance of coherence (Landau- 
Pomeranchuk [TT] and Migdal \12 ) effects for large enough 
formation times of the emitted gluon, so that Tform^Amfp- 

If the soft-gluon radiation were the only mechanism re- 
sponsible for energy loss, one would expect heavy quarks 
(c and b) to be much less quenched. Characteristic ef- 
fects occurring for massive particles are the suppression 
of coUinear gluon emission for angles 9<M/E {dead cone 
effect [I31II1]), the shortening of the gluon formation time 
and the suppression of the radiated gluon spectra at large 
energy [15 . Furthermore, while light-hadron spectra get 
contribution from the fragmentation of partons belonging 
both to the fundamental (quarks) and the adjoint (gluons) 
representations, heavy-flavour electrons arise only from 
the decays of c and b quarks. Hence, the associated Casimir 
factor Cp (smaller then Ca of the gluons) should entail 
a larger mean-free-path, resulting in a lower rate of gluon 
radiation. However, single-electron spectra (resulting from 
the semi-leptonic decays of D and B mesons) measured by 
the PHENIX experiment at RHIC show a similar level of 
suppression as the one found for light hadrons [TBlfTT] (a 
similar analysis is being performed by the STAR collabo- 
ration); moreover, they also display a sizable elliptic flow. 
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These findings hint at reconsidering the importance of col- 
hsional energy loss. Calculations of the stopping power of 
the quark-gluon plasma (QGP) for heavy quarks due to 
elastic collisions can be found for instance in Refs. [TFIITOI 

An appropriate tool to study the final spectra of heavy 
quarks, after their evolution in the fireball created in the 
heavy ion collisions, is given by the (relativistic) Langevin 
equation, which relies on the picture of many uncorrelated 
random collisions |21[|22| . This approach — pursued by 
solving either the Langevin equation or the correspond- 
ing Fokker-Planck equation — has already been explored 
in the literature P51IM1[^[^[?71[^ and will be the one 
addressed in this paper. Solving the Langevin equation re- 
quires as an input the knowledge of a few transport coeffi- 
cients, accounting for the interaction of the heavy quarks 
with the medium: indeed, in this framework, the latter 
gives rise to a friction force and a momentum broaden- 
ing. The various calculations available in the literature 
[11IM1I251I2S1I2Z1I1H] differ mainly in the way they esti- 
mate the above coefficients, either inspired by pQCD or 
within less conventional scenarios (assuming for instance 
the existence of resonant states in the plasma |23] or the 
AdS/CFT correspondence In all the cases the free- 

dom of tuning some parameter allows one to explore more/l( 
strongly coupled regimes. 

However, before looking for alternative explanations, 
it appears important to have at hand a fully consistent 
weakly coupled calculation, to check how closely the ex- 
perimental data can be reproduced within a perturbative 
picture and how much room is left for more "exotic" sce- 
narios. This is the major goal of our paper, in which we 
provide a microscopic evaluation of the heavy-quark trans- 
port coefficients and we use them to solve the Langevin 
equation for c and b quarks that are produced in nucleus- 
nucleus collisions and are let propagate randomly in the 
resulting deconfined medium. Our calculation relies on the 
separate treatment of soft and hard collisions. The first 
ones, mediated by the exchange of long wavelength gluons, 
will be described within the Hard Thermal Loop (HTL) 
approximation, which resums medium effects. Hard colli- 
sions, involving a high-momentum exchange, will be cal- 
culated in kinetic theory, employing leading order pQCD 
matrix elements The two contributions will then be 
summed up. Once the transport coefficients have been 
calculated, we can address the numerical solution of the 
Langevin equation, given suitable initial conditions and a 
background medium. 

The initial heavy-quark pairs are produced using the 
POWHEG package [501I5T] . a hard event generator for 
heavy-quark production in hadronic collisions, which im- 
plements pQCD at next-to- leading order (NLO). Since it 
generates events, this code is particularly suitable for a 
Langevin simulation. Care is taken to allow also for shad- 
owing and transverse momentum broadening effects. Once 
created, the heavy quark is allowed to propagate in the 
hot medium according to the Langevin dynamics, with 
transport coefficients calculated as mentioned above. The 
evolution of the background medium is described by two 



Table 1. Total production cross sections for heavy-quark pairs 
calculated with POWHEG for pp collisions at two different 
center-of-mass energies and for two choices of PDF's; CTEQ6M 
alone (suitable for bare pp collisions) and supplemented by the 
EPS09 nuclear modifications 39 (for NN events embedded in 
a nucleus-nucleus collision). 





(Tec (mb) per NN event 


(T^j, (mb) per NN event 




CTEQ6M (+EPS09) 


CTEQ6M (+EPS09) 


200 GeV 
5.5 TeV 


0.254 0.236 
3.015 2.288 


L77xl0~^ 2.03x10"^ 
0.187 0.169 



different scenarios, namely ideal [5^1551151] and viscous 
35,36,37 hydrodynamics, in order to estimate the de- 
pendence of our results on the hydrodynamical scheme. 
Around the phase transition energy density the c and b 
quarks are made hadronize (at the moment, for the sake 
of simplicity, only through fragmentation) and then decay 
to electrons, to allow for a comparison with the existing 
experimental data from RHIC for open heavy-flavour elec- 
trons. Results will be presented for the invariant yields, 
the nuclear modification factors and the elliptic flow coef- 
flcients at the kinematics of both RHIC and LHC. For the 
;Jatter we have chosen the highest center-of-mass energy 
that the heavy ion experiments should be able to reach, 
i. e. ^/s = 5.5 TeV. Preliminary results have already been 
presented in Refs. [55] . 

Our paper is organized as follows: in Sect. [5] we in- 
troduce the theoretical tools employed in generating the 
heavy-quark spectra after a high-energy heavy ion col- 
lision, examining all the steps of the process, from the 
production stage, through the hydrodynamic evolution, 
the Langevin propagation, the calculation of the transport 
coefficients, up to hadronization and decay; in Sect. [3] we 
present and discuss our findings for the differential spec- 
tra, the nuclear modification factor and the elliptic-fiow 
coefficient; finally, in Sect. H] we sum up our conclusions; 
some more technical material is deferred to the Appen- 
dices. 



2 Theoretical framework 

2.1 Heavy-quark production in pp and AA collisions 

For every experimental setup (either pp or AA at ^/s — 
200 GeV (RHIC) and = 5.5 TeV (LHC)) we have gen- 
erated 4.5 X 10^ heavy-quark pairs (either cc or bb) using 
the POWHEG event generator [5011311 . which implements 
pQCD at NLO. 

There is considerable uncertainty in the QQ produc- 
tion cross section, due to uncertainties in the quark masses 
and in the renormalization and factorization scales. Here 
we have chosen rric = 1.5 GeV, rrib = 4.8 GeV and for 
the values of the renormalization and factorization scales 
we have kept the default choice in POWHEG, that is the 
heavy-quark transverse mass in the QQ reference frame 
[50] . For pp events we have employed the CTEQ6M parton 
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Fig. 1. (a) Invariant differential cross section of electrons 
[(e+ -I- e~)/2], from heavy-flavour decay in pp collisions at 
= 200 GeV, at mid-rapidity as a function of transverse mo- 
mentum. The experimental points represent PHENIX data [421 
I17| ; both statistical (error bars) and systematic (gray boxes) er- 
rors are displayed. The curves represent the contributions com- 
ing from electrons originating from a charm or bottom quark 
generated by POWHEG with inclusion of transverse broaden- 
ing, (b) ratio of data to the POWEIEG prediction as a function 
of pt, with (circle) and without (square) inclusion of transverse 
broadening; the lines represent the best fit to a constant. 



distribution function (PDF). The total cc and bb produc- 
tion cross sections one gets at RHIC energies (see Table [1]) 
are very close to the central values predicted by FONLL 
[40] . Another important effect to be accounted for in the 
production cross sections is related to the intrinsic trans- 
verse momentum of the initial partons. We have adopted 
the procedure described in Ref. [JT] , which amounts to add 
to the out-coming heavy quarks a transverse momentum 
contribution randomly generated from a Gaussian distri- 
bution with given variance (A;|,)nn- For pp collisions we 



have adopted the standard value (fc^) 



NN 



1 GeVVc^ 



In Fig. [T] one can see the invariant differential cross sec- 
tions of electrons from heavy-flavour decay (see Sect. 12.51 
for details about hadronization and decay) in pp collisions 
at y/s — 200 GeV, compared to the mid-rapidity data 
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Fig. 2. Invariant differential cross section of electrons [(e^ -|- 
e~)/2], from heavy-flavour decay in pp collisions at y^s = 
5.5 TeV, at mid-rapidity as a function of transverse momen- 
tum. The curves represent the contributions coming from elec- 
trons originating from a charm or bottom quark generated by 
POWHEG with inclusion of transverse broadening. For com- 
parison we have also included the total yield of electrons at 
RHIC from Fig. [T] 



from PHENIX [3^[T7pl . One observes that the pQCD out- 
come provides a fairly good description of both the shape 
and the absolute magnitude of the data. In panel (b) of 
the figure the ratio between experiment and theory is also 
displayed: this ratio is nearly pT-independent and a best 
fit with a constant provides for the latter the value 1.12. 
Given a 10% normalization uncertainty in the data (not 
shown in the figure) and the previously mentioned uncer- 
tainties on the theoretical parameters, a 12% difference 
between theory and data can be easily accommodated. Ac- 
cordingly, for sake of simplicity, in the following we have 
chosen to multiply by this factor all the calculated cross 
sections at RHIC energies (of course, Raa and V2 are not 
affected by cross section normalization). In Fig. [Ud we 
also display the ratio of data to the cross section gen- 
erated without any transverse momentum broadening: in 
this case the discrepancy is fairly larger (although still 
within the experimental and theoretical uncertainties). 

In Fig.[2]the invariant differential cross sections of elec- 
trons from heavy-flavour decay in pp collisions at \fs = 
5.5 TeV are reported: apart from the much larger strength 
and the less steep slope of the cross sections, one observes 



^ For the purpose of comparing to the PHENIX data, the 
doubly differential cross sections have been calculated by in- 
tegrating over the experimental acceptance in pseudo-rapidity 
(|?7| < 0.35) and dividing by the rapidity range (Z\y = 0.7). For 
consistency, also the cross sections at the energy of LHC have 
been calculated in the same way, using the constraints of the 
ALICE experiment (|r;| < 0.9 and Ay ~ 1.8). 
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that the bottom relative contribution is much more im- 
portant than at the RHIC energy. 

In the case of nucleus-nucleus collisions, there are two 
important differences one has to consider in connection to 
the initial heavy-quark distributions. First of all, the nu- 
clear PDF's should be different from the ones employed 
in pp collisions and, to account for this occurrence, we 
have adopted here the EPS09 scheme [39j. In principle, 
the density probed by the colliding partons should depend 
on the impact parameter b: in describing nucleus-nucleus 
collisions we have made the simple choice of employing 
the EPS09 scheme for impact parameters b < 2R and of 
neglecting nuclear corrections for b > 2R {R being the 
radius of the nuclear density distribution) . The main con- 
sequence of using a different PDF in AA and pp collisions 
relates to the different total (per binary collision) cc and 
bb production cross sections one gets, as one can see in 
Table [TJ Indeed, since in a binary process the longitudi- 
nal momentum fraction carried by the two initial partons 
(most of the times gluons) is given by [13] 



3^1/2 = MqQ 



NN 



(1) 



Mqq and Hqq being the invariant mass and rapidity of the 
QQ pair, one can get larger or smaller cross sections in A A 
collisions, at different energies and quark masses, depend- 
ing upon whether one is probing the anti-shadowing or 
shadowing regions. For instance, at RHIC energy, because 
of the large bottom mass, the main contributions come 
from relatively large values of x, where anti-shadowing 
dominates — hence a larger bb cross section per binary 
collision results — whereas at LHC energies one is mainly 
probing the low x shadowing region, yielding smaller cross 
sections. This fact has sizable consequences on Raa, as we 
shall see later. 

A second difference one has to cope with in AA colli- 
sions concerns the larger transverse momentum acquired 
on average by the colliding partons, because of the large 
size of the traversed medium. To get a realistic estimate 
for {k'^)AA in nucleus-nucleus collisions we have adopted 
a Glauber approach — developed in Refs. in or- 

der to study the pr distribution of charmonia produced in 
proton-nucleus and nucleus-nucleus collisions [46] — and 
we have generalized it to study the inclusive single-quark 
spectra. Details of the procedure can be found in Ap- 
pendix|21 One gets an average squared transverse momen- 
tum that depends not only on the impact parameter of the 
collision and on the nuclei involved, but also on the trans- 
verse position of the QQ pair. To get an orientation on the 
magnitude of the quantities involved, we mention that at 
RHIC energy one gets, depending upon the impact param- 
eter, average values of {k^)AA around 1.3 1.5 GeV^/c^ 
for charm and around 1.8 2.3 GeV^/c^ for bottom in 
Au-Au collisions; at the LHC {y/s — 5.5 TeV) around 
1.5-^1.8 GeVVc^ for charm and around 2.5^3.3 GeVVc^ 
for bottom in Pb-Pb collisions. 

Finally, for A A collisions, since we are following in 
the Langevin framework the space-time propagation of an 
ensemble of heavy quarks, we have also to speciiy their 
space-time distribution at the onset of the hydrodynami- 



cal evolution of the background medium. This we do con- 
sistently with the choice discussed below for the initial 
conditions of the hydrodynamic equations. Namely, at the 
initial proper time tq the position of the heavy quarks in 
the transverse plane is calculated in a Glauber framework 
according to a distribution generated by the nuclear over- 
lap function T{x + 6/2, y)T{x ~ b/2, y), where 



Tix,y) = / dzp{x,y,z), 



(2) 



p being a Fermi parameterization of the nuclear density 
[47) and b the impact parameter. In the longitudinal di- 
rection we set z — Tq sinh r]s , with 

,..ilni±i = iln^. (3) 
' 2 t-z 2 E-P;, ^ ' 



2.2 Hydrodynamic evolution at RHIC and LHC 

Hydrodynamics has been successfully applied to the de- 
scription of collective phenomena in heavy ion collisions 
at RHIC, yielding a sensible description for a number of 
experimental observables. For our purposes, hydrodynam- 
ics provides the full space-time evolution of the properties 
of the expanding medium — such as temperature, flow 
velocity and energy density — that are needed to follow 
the propagation of the heavy quarks. We have chosen two 
different implementations of the relativistic hydrodynamic 
equations, whose codes are publicly available, namely ideal 
[5^[551[51] and viscous [5^1351137] hydrodynamics, both as- 
suming exact longitudinal boost invariance. The two mod- 
els differ not only in the ideal/viscous implementation 
(with the ratio of shear viscosity to entropy density taken 
to be rj/s = 0.08 in the viscous case), but also in the 
choice of the equation of state (EOS) and of the initial 
conditions. By comparing the results obtained in the two 
scenarios one can get an estimate of the amount of uncer- 
tainty due to the treatment of the background medium. 

The initial energy density distribution is computed in 
both cases according to the Glauber model, using either 
the number of participating nucleons (A^part) or the num- 
ber of binary collisions (A'coii)- In the ideal model [S^ 
the distribution has been ascribed for the 75% to "soft" 



Table 2. Values of the initial proper time to, central entropy 
density sq and central temperature To employed as initial con- 
ditions for the hydrodynamical evolution. 





■q/s = 


ri/s = 0.08 




TO (fm) So (fm--')ro (MeV) 


TO (fm) so (fm-')To (MeV) 


RHIC 


0.6 110 357 


0.1 840 666 
0.6 140 387 
1 84 333 


LHC 


0.1 2438 1000 
0.45 271 482 


0.1 1840 854 
1 184 420 
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processes scaling as iVpart and for the 25% to "hard" pro- 
cesses scaUng as A^coii; hi the viscous model the initial en- 
ergy density is proportional to iVcoii- Concerning the EOS, 
in the ideal case it is obtained by matching, through a 
Maxwell construction, a gas of hadron resonances with an 
ideal gas of massless quarks and gluons; this model yields 
a first order transition at Tciit = 164 MeV. In the vis- 
cous case, one employs the more realistic EOS of Ref. ^48j, 
which implements a crossover between a low temperature 
hadron-resonance gas and a high temperature pQCD cal- 
culation. 

The initial proper time tq is one of the parameters of 
the hydrodynamical model: at the conditions of RHIC, 
the values tq = 0.6 fm and tq = 1 fm have been chosen 
by the authors in the ideal [32] and viscous [5B| cases, re- 
spectively. Since the value of tq has some impact on the 
final heavy quark spectra (a lower value provides a longer 
propagation and exposes the heavy quark to higher tem- 
peratures), we have also explored a scenario with a very 
small value of tq, by adjusting the maximal initial entropy 
density Sq (hence, through the EOS, the corresponding ini- 
tial energy density) in such a way to maintain the same 
particle density, using the relation 



dN 
dy 



(4) 



At the LHC regime, on the other hand, a full experimen- 
tal constraint is not available yet (although first data on 
particle density at mid-rapidity are becoming available at 
= 2.76 TeV 09]). Here we have followed Refs. [50l 
[CT] (ideal) and Ref. [37| (viscous), where the initial con- 
ditions have been fixed in order to match, instead of the 
experimentally observed multiplicities, the range of the 
predicted ones. Our set of initial conditions is summed up 
in Tabled 

As an example of the differences among the various hy- 
drodynamical scenarios that we have employed, in Fig. |3] 
we display the temperature and the energy density at the 
center of the fireball during the hydrodynamical evolution 
for a few cases of semi-peripheral collisions. For every ex- 
perimental setup, at RHIC and LHC, we have reported 
the scenarios giving rise to the lowest and to the highest 
temperatures attained during the evolution. 



2.3 Langevin dynamics in a relativistic fluid 

We wish to follow the propagation of c and b quarks — 
initially produced in hard pQCD processes — in the ex- 
panding, thermalized fireball of deconfined matter that 
one expects to arise from a high energy heavy ion collision. 
The propagation of the heavy quarks in such a hot envi- 
ronment is modeled as a Brownian motion by employing 
a relativistic Langevin equation [21,23,24,26 : our imple- 
mentation of the model has been discussed in detail in a 
previous paper [22], where we dealt with the simpler case 
of a static homogeneous medium at rest. Here, we briefly 
summarize the essential points of model. 




Fig. 3. (a) The temperature at the center of the flreball dur- 
ing hydrodynamical evolution for a few scenarios both at RHIC 
(impact parameter h = 8.44 fm) and LHC {b = 8.77 fm). (b) 
The energy density at the center of the fireball during hydro- 
dynamical evolution for the same scenarios. 



The evolution with time of the momentum of a rel- 
ativistic Brownian particle is provided by the Langevin 
equation: 

dp 



-± = -vDiP)p+m, 



(5) 



where the drag coefficient rjoip) describes the determin- 
istic friction force acting on the heavy quark, whereas 
the term ^ accounts for the random collisions with the 
constituents of the medium. The effect of the stochastic 
term is completely determined once its temporal correla- 
tion function is fixed. The latter is usually assumed to be 
given by 



{C{tW{t'))^V^{p)5{t-t')^ 



(6) 



entailing that collisions at different time-steps are uncor- 
related. The tensor (p) can be decomposed with a stan- 
dard procedure according to 



(7) 



in terms of the coefficients KL/rip), which represent the 
squared longitudinal/transverse momentum per unit time 
exchanged by the quark with the medium. It is useful to 
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introduce the related tensor 



2.4 Heavy-quark transport coefficients 



^9L{p)fi^ +9t{pW -ff), (8) 

which ahows one to factor out the momentum dependence 
of the noise term in Eq. (O , thus obtaining the equation 



= -vd(p)p' +g'\p)rf{t), 



(9) 



with 



Wit)vHt'))^s^'S{t~t'). (10) 

Finally, the drag coefficient rju^p) is fixed in order to en- 
sure the approach to equilibrium: for large times the mo- 
menta of an ensemble of heavy quarks should be described 
by an equilibrium Maxwell- Jiittner distribution. Actually, 
the whole procedure depends on the discretization scheme 
employed in the numerical solution of Eq. (|9l), which be- 
longs to the class of the stochastic differential equations 
[5^155] . In the Ito discretization scheme one gets: 



i^l{p) 
2TE 



1 



(1- 



Qy2 y2 



(11) 

where v is the quark velocity, E its energy and T the 
medium temperature. 

The set of Eqs. (|51 fTT|) is defined in the rest frame 
of the background medium and it allows one to study 
the quark propagation once the transport coefficients are 
given. These depend on the medium temperature, which in 
turn, in the expanding fireball, depends on the space-time 
position occupied by the heavy quark. Hence, to study the 
fate of a heavy quark through the quark-gluon plasma, we 
adopt the following procedure: 

— We determine the initial four-momentum p^^ and the 
initial space-time position a;'' of the heavy quark (in 
the laboratory system) as explained in Sect. 12. ll 

— Given the position x'^, we use the information from 
the hydrodynamic simulation (Sect. [^3) to obtain the 
fluid local temperature T{x), velocity u^{x) and en- 
ergy density e{x). 

— We check (see Sect. 12. 5p whether the conditions for 
hadronization apply: in this case the procedure is ended; 
otherwise 

— we make a Lorentz transformation (p^ p^) to the 
fluid rest frame, employ Eqs. (|51 ITT|) to update the 
quark momentum (p^ p'^) and boost it back to 
the laboratory (p''' — > p'^). 

— We update the space-time step made by the quark in 
the rest frame {Axf^ = {p'^/Ep)Ai), boost it to the 
laboratory {Ax^ — >■ Ax^) and use it to update the 
quark position (a;'' a;'^). 

— Given the new momentum p''^ and the new position 
x''^ the procedure is started again until the conditions 
for hadronization are met. 

The time step in the rest frame, which enters in updating 
the quark position and also the quark momentum through 
the Langevin equation, in our calculations has the value 
Z\t = 0.02 fm. 



In the approach to stochastic dynamics based on the (rela- 
tivistic) Langevin equation, the information on the prop- 
erties of the medium crossed by a Brownian particle is 
encoded in the drag and momentum-diffusion coefficients 
r]D and k-t/l- Following the random motion of the heavy 
quarks in the QGP requires then the calculation of the 
above coefficients, starting from the microscopic theory, 
i.e. finite temperature QCD. In accord with Refs. P^I20) . 
we accomplish this task by separating soft and hard col- 
lisions depending upon the value of the Mandelstam vari- 



able t 



q that characterizes the elementary scat- 



tering, and summing at the end the two contribution^ 

Soft scatterings — corresponding to i/jij ~ gT — oc- 
cur quite frequently, with a typical mean free path X^^fp ~ 
1/g^T. Being related to the exchange of long- wavelength 
gluons, they require a careful treatment of medium effects, 
which can be achieved by describing these soft collisions 
within the HTL approximation. Hard scatterings — corre- 
sponding to \/\t\^T — are more rare, the mean free path 
being in this case X^lp ^ l/g'^T, but, since they cause 
a sizable momentum exchange, their contribution to the 
transport coefficients is substantial and will be evaluated, 
as already mentioned, within a microscopic kinetic calcu- 
lation based on pQCD. We set \t\* ^ m|), niD being the 
Debye mass, as the intermediate cutoff. Clearly, the final 
results should not be too much affected by the choice of 
\t\* and we have verified that this is actually the case, 
even for temperatures of experimental interest, where the 
coupling is not really small. 

The coupling constant g has been chosen to run ac- 
cording to the two-loop QCD beta-function, with vIqcd — 
261 MeV, as in Ref. [S3]. The coupling constant displays 
a strong dependence on the unknown value of the tem- 
perature dependent scale fi (x T and this has important 
consequences on the final observables, as we shall see in 
the following. To shed light on this important issue, we 
shall explore the effect on the calculated quantities for a 
wide range of values, ttT < /i < 2ttT. 



2.4.1 Soft collisions: hard thermal loops 

The contribution of soft collisions to the transport coeffi- 
cients has already been obtained in our previous work [35] , 
to which we refer for details. However, due to our choice 
of setting \t\* as an intermediate cutoff, it is necessary 
to rewrite the formulas in Ref. \22\ for Kt/l employing 
\t\ = cj^ — uP' and x = ui/q as integration variables. One 
gets: 



.soft 



tl* 



d\t\ 



dx—— 
2(1 



2)5/: 



rPd^U) 



Note that this differs from what is usually done in the liter- 
ature (see, e.g., Ref. [H]), where the separation between hard 
and soft scatterings is related to the spatial momentum q ex- 
changed in the collision. 
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Fig. 4. The tree level diagrams for the hard scattering of a heavy quark off a light (anti-)quark of the thermal bath and for the 
hard scattering (in the t, s and u channels) of a heavy quark off a gluon. 



xll-^ cothp^^L 1\ (12) 



and 



^soft 



Att^v Jo ' ' 7o 2(1 - ^2)5/2 



pi\t\,x) 



I x^\t\/{\--x-^) 



X — coth 



2T 



(13) 



for the transverse and longitudinal momentum diffusion 
coefficients, respectively. In the above, v is the quark ve- 
locity, Cp — 4/3 the quark Casimir factor and 

p{\t\,x) = PL{\t\,x) + (v^ - X^)pT{\t\,x), (14) 

where Pl/t ^^re the continuum parts of the HTL gluon 
spectral functions: 



PL{\t\,x) = Trmjyx 



l-x^ 



2 / 1 *^ 1 



PT{\t\,x) = 7rm|,-(l - x^) 



X + 


1 




2 

+ 




x — 


1 


)] 





1 - j;2 , 



x+l 



x-l 



\t\+mi,- 



" 2-^/1 2 
— mn— 1 — X 

2 ^2^ 



(15) 
\ -1 



2.4.2 Hard collisions: perturbative QCD 

The contribution from hard collisions to the transport co- 
efficients ht/l can be calculated starting from their mi- 
croscopic definition, namely 



1 / Aq^ 



and kl 



At 



(16) 



and by weighting the interaction rate with the squared 
transverse and longitudinal momentum exchanged in the 
collisions with gluons and (anti-)quarks of the medium. 
One has 



where (employing the notation J^, = J dk/{27T) 



/g_ I f riB/pik) f 1 ± riB/Fik') f 9i\t\ ~ \t\*) 



2E ./fe 2k 



2k' 



p' 



2E' 



X (2^)4^(4) ^p^K-P'- K') \M,,, {s, r ^ (18) 




10 15 
p (GeV/c) 



20 



Fig. 5. The charm (upper panel) and bottom (lower panel) mo- 
mentum diffusion coefflcients kt/l{p) resulting from the sum 
of the soft and hard contributions. The sensitivity to the in- 
termediate cutoff ~ mjj is very mild. The curves refer to 
the temperature T = 400 MeV, with the coupling g evaluated 
at the scale p = (3/2)7rT. 
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and 



1 /■ nB/F{k) 



2k 



l±nB/F{k') f e{\t\-\t\*) 



2k' 



p' 



2E' 



X (2^)4^(4) (P + K-P'-K') \M,,,[s, t)\^ ql (19) 

In the above expressions, s = (P + K)^, Q = (w, q) = P — 
P' 1 1tI<Il are the transverse/longitudinal components of 
q and n^/F are the Bose/Fermi distributions; the squared 
ampHtudes are averaged over the internal degrees of free- 
dom (colour and spin) of the incoming heavy quark and 
summed over the internal degrees of freedom of the ini- 
tial and final state partons in the medium: A^^ — 1 colours 
and two transverse polarizations for a gluon, TV^ colours 
and two helicities for a light quark. Furthermore a fac- 
tor 2N f is assigned to , accounting for the identical 
contribution provided by all the flavours of light quarks 
and anti-quarks in the medium. The squared amplitudes 

|A^g/g| defined in such a way can be obtained starting 
from the results given in Ref. [29] and their expressions 
are reported in Appendix [BJ The quark contribution is 
then obtained by squaring the corresponding amplitude in 
Fig. |4l given by a simple gluon exchange in the i-channel. 
The gluon contribution is more cumbersome and requires 
squaring the sum of the s and u-channel diagrams in 
Fig.H 



2.4.3 Charm and bottom transport coefficients 

In this section we display the outcome of our calculation 
of the heavy quark transport coefficients in the QGP and 
compare our results with previous findings for the same 
coefficients available in the literature [22.231 IM1I25) . 

As we have already mentioned, following the heavy- 
quark stochastic propagation in the hot environment pro- 
duced in heavy ion collisions requires the knowledge of the 
corresponding transport coefficients for a range of tem- 
peratures from T ^ Tc (where actually we expect our 
weak-coupling approach to underestimate their value) up 
to T ~ 1 GeV (for the case of the LHC scenarios with the 
highest initial energy-density). 

In Fig. [5] we show our findings for kt/l(p) for c and 
b quarks, at the temperature T — 400 MeV; the coupling 
g has been evaluated at the scale /i = (3/2)7rT. Two im- 
portant features are apparent in Fig. [5] the larger growth 
with the momentum p of kl{p) with respect to ht{p) and 
the very mild sensitivity to the value of the intermediate 
cutoff ~ ml,. This last finding occurs in spite of the 
fact that, at the experimentally accessible temperatures, 
the coupling is not small. Such an occurrence supports the 
validity of the adopted approach. 

A deeper insight can be gained by looking at Fig. [51 
where we plot separately the soft and hard contributions 
to Krip) and kl{p) for the case of c quarks. 

We finally compare our results for the heavy-quark 
transport coefficients to the ones obtained using other 
models. We start with our previous work [22], where the 
HTL approximation had been applied to any value of the 




p (GeV/c) 

Fig. 6. The separate soft and hard contributions to ht{p) 
(upper panel) and kl{p) (lower panel) for different choices of 
the intermediate cutoff \t\* ~ mf,. The curves refer to the case 
of a c quark, in a medium at T = 400 MeV, with the coupling 
g evaluated at the scale /i = {3/2)nT. 



momentum exchange, up to a cut-off (/max- Results for c 
and b quarks are shown in Fig. [7^ and compared with the 
present, more realistic, treatment. The major difference 
one observes is in the growth of kt{p) with p, which is 
faster in the pure HTL approach. 

In Fig. [TJa we summarize the outcomes of several mod- 
els which can be found in the literature: the present ap- 
proach (pQCD+HTL), the pure HTL result [22] and the 
findings of Refs. 23,24,25j. In the approach of Ref. |23) . 
the coefficients are obtained from a kinetic pQCD calcu- 
lation, with the value of k{p = 0) tuned by hand in order 
to explore both strongly and weakly-coupled scenarios. In 
Ref. 1 23^ the mechanism responsible for the heavy-quark 
thermalization is resonant scattering, with the temporary 
formation of finite- width D-mesons in the plasma. The au- 
thors of Ref. [53] J on the other hand, resort to an effective 
running coupling constant and to a value of the Debye 
mass lower than the thermal QCD prediction. Notice how 
all the models based on pQCD tend to favor small-angle 
scattering. 
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Fig. T. (a) Transport coeflicients, for both c and b quarks 
at T = 400 MeV and fj. = {3/2)ttT, arising from the sepa- 
rate treatment of soft (HTL) and hard (pQCD) collisions (with 
= m|,) compared with the pure HTL results of Ref. 22:. 
(b) Comparison of our results for c quarks (with \t*\ = m|, 
and jj, = {3/2)ttT) to the outcomes of several models [2211231 
[241125] for the transport coefficients kt/l{p) at T = 200 MeV. 



2.5 Hadronization and decay 

The conversion of the heavy quarks into hadrons requires 
two distinct steps: first, one has to decide when the quark 
is going to hadronize and then to apply a specific model 
for the transition. In the ideal hydrodynamics model of 
Refs. a first order phase transition occurs at 

the temperature Tc = 164 MeV, lasting for a couple of fm 
and during which the energy density drops from egcp = 
1.65 GeV/fm3 to Eh = 0.45 GeV/fm^. In the viscous hy- 
drodynamics model of Refs. [3511361137] , on the other hand, 
a more realistic crossover, slightly faster because of viscos- 
ity, is employed. In Fig. [S] we display the energy density in 



1.5 - 



I Ih 
> 

o 



0.5 



J I I y I i\ I I I I I I I y I r^i I I I I I I I I I I I I I I I 



■\ ■ n 

\ \ 
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\ 



\ 

\ 

\ 



RHIC, ideal, x=0.6 fm 

RHIC, viscous, T=l fm 

LHC, viscous, x=0.1 fm 

LHC, viscous, x=l fm 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



3 4 5 6 7 8 9 
X(fm) 

Fig. 8. Energy density in the center of the fireball at proper 
times around the deconfinement transition for a few hydro- 
dynamical setups; the dotted lines correspond to eqgp = 
1.65 GeV/fm^ and eh = 0.45 GeV/fm^ 



the center of the fireball around the QGP/hadron- matter 
transition for a few hydrodynamical scenarios. The dotted 
lines mark the position of £qgp and en: we assume that 
hadronization is going to take place in this region of mixed 
phase. 

Introducing the fraction of QGP in the mixed phase 
pS] according to 



/qgp 



en 



£qgp — Eh 



(20) 



we stop the Langevin propagation of the heavy quark ac- 
cording to the following prescription: 

— we extract the medium energy density at the heavy- 
quark space-time position; 

— if /qgp is larger than one, the Langevin propagation 
is carried on another step; otherwise 

— we treat 1 — /qgp as a transition probability: given 
/qgp — fraction of QGP at the i-th propaga- 
tion step — we draw a random number h; if /i > 
/qgp//qgp ^'^^^ hadronization has occurred, other- 
wise a new Langevin step is taken and the procedure 
repeated at the new position. 

In this way the transition from quarks to hadrons is made 
occur over the whole mixed phase. 



2.5.1 Hadronization 

Two mechanisms of hadronization should in principle be 
implemented, namely recombination of the heavy quark 
with a light quark from the medium and fragmentation. 
The coalescence approach has been shown to be important 
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2 4 6 8 02468 10 12 



(GeV/c) (GeV/c) 

Fig. 9. Nuclear modiflcation factor Raa vs. pr for open charm (left) and open beauty (right) hadrons for different values of 
the e parameter in the Peterson fragmentation function; the dotted lines represent the delta function limit in the fragmentation 
function, D{z) = 5(1 — z); in these examples we have used ideal hydrodynamics at RHIC for b = 8.44 fm and /i — 2nT. 



in the low momentum part of the spectra at RHIC 
say for pt^3 GeV/c. In this work, we limit ourselves to let 
the heavy quarks hadronize by fragmentation, hence our 
results should be viewed as realistic for relatively large 
momenta, ptJ^3 GeV/c at RHIC energy. Note, however, 
that the region where recombination is dominant might 
change at LHC energies. 

Open charm and open beauty hadrons are produced 
from the charm and bottom quarks via a two-step Monte 
Carlo procedure. First, the resulting hadron species is as- 
signed according to measured fragmentation ratios. In the 
case of charm quarks, four hadronic final states have been 
considered, namely D°, D^, D+ and A^, together with 
their corresponding anti-particles. The fragmentation ra- 
tios are taken from Table 4 of Ref. [SS] . For bottom quarks, 
the considered hadronic states have been B", B+, and 

together with their corresponding anti-particles. The 
fragmentation ratios have been taken from with the 
only modification of assigning to the yl" final state the 
probability for a bottom quark to fragment into b-baryons. 

The second step is the Monte Carlo generation of the 
hadron momenta. This has been done starting from the 
momentum components (px,Py and Pz) of the heavy quark 
out-coming from the Langevin evolution and sampling, 
from a Peterson fragmentation function [57] , the fraction z 
of quark momentum that is taken by the produced hadron: 

D{z)=^--(l-^---^Y . (21) 
z \ z 1 — z J 

The parameter k has been set in order to have the frag- 
mentation function normalized to 1, while for e the val- 
ues 0.04 and 0.005 have been used for charm and bottom 
quarks, respectively. They have been chosen by compar- 
ing the fragmentation function shapes with the ones com- 
puted in Ref. [55] on a pQCD or Heavy Quark Effective 
Theory basis. To check the effect of the choice of the e pa- 
rameter on the resulting nuclear modification factor Raa 
(see Ea. E5|) . the whole fragmentation procedure has been 



repeated four times varying the value of e in the range 
yielding a reasonable agreement with the calculation of 
Ref. [52 . The obtained Raa for open charm and open 
beauty hadrons as a function of the transverse momen- 
tum pt for different choices of e is shown in Fig. [5] for 
a specific case of -^/s, centrality and hydrodynamical sce- 
nario. We have checked that variations of the parameter 
£ by a factor two results generally in a <3% effect on the 
Raa of open heavy-fiavoured hadrons. 

2.5.2 Decays of heavy-flavour hadrons into electrons 

In order to compare our calculations to the experimental 
data for the non-photonic single-electron (e^) transverse 
spectra and for the corresponding nuclear modification 
factor Raa and elliptic fiow coefficient W2, each charm or 
bottom hadron originated by the fragmentation of a heavy 
quark (as described in the previous section), is forced to 
decay into electrons. 

This task is performed by using the decay routines im- 
plemented inside the P YTHIA event generator [SU] • How- 
ever, a revised table of the branching ratios for the domi- 
nant decay channels of each heavy-quark hadron has been 
employed, in order to replace the default values which are 
by now obsolete. Such table has been drawn up by taking 
the measured values of the branching ratios, as reported 
on the latest Particle Data Group review [50] . 

Only hadronic and semi-leptonic decay channels of charm 
and beauty hadrons with branching ratio F > 10^^ have 
been included. Therefore the contribution of each hadron 
to the final electron spectra has been properly weighted 
by the corresponding total branching ratio. 

In this way, all contributions to the electron yield (i. e., 
prompt electrons from i? — > e or I? — >■ e decays, and 
secondary electrons from a cascade process i? — > 13 — )■ 
e) are naturally implemented in our simulation scheme, 
as well as the additional "background" contribution from 
the decay chain B J/t/j e~^e~. Indeed, the latter 
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contribution has been subtracted from the total measured 
electron yields only in the most recent analyses of the 
PHENIX data [TT]. 

The electron spectra originated from charm (cc) and 
bottom (ch) decays are then combined into a unique spec- 
trum (cc+b) with an appropriate weight, which accounts 
for the corresponding production cross section of the par- 
ent heavy quark. 



3 Heavy-flavour spectra in nucleus-nucleus 
collisions 

We have considered two kinematic regimes for high energy 
heavy ion collisions: Au-Au at y/s = 200 GeV as in the 
PHENIX experiment at RHIC dldT] and Pb-Pb at = 
5.5 TeV, the highest energy that should be attained with 
A A collisions at the LHC. 

In the PHENIX experiment the heavy-flavour data 
have been presented in several centrality classes C1-C2, 
corresponding to the fraction of the geometric cross sec- 
tion with impact parameter in the range 61 < b < b2- 



fCi-C2 



JI^^ dbb[l - exp(crNN7'yis(&)] 
dbb[l - exp(crNNTAB(fe)] 



(22) 



where Tab(6) = / dsTA{s + b/2)TBis-'b/2) is the nuclear 
overlap function and the nuclear profile function in the 
transverse plane Ta{s) is given in Eq. ([2]). Our calculations 
are performed for a given impact parameter 6, so for every 
centrality class we have to evaluate an average impact pa- 
rameter. This can be done by equating the number of col- 
lisions in a given centrality class < iVcoii >Ci-C2— "'nn < 
Tab >Ci-C2= o'nn JI^^ dbbTAsib)/ dbb to the num- 
ber of collisions at a given impact parameter iVcon(6) — 
i^NNTABib) [61]. In Table [3] we show the centrality classes 
and the corresponding impact parameters that we have 
considered in the present work. At the LHC energy, lack- 
ing the experimental data, we have kept a partition similar 
to the one employed of RHIC. 



3.1 Differential spectra 

In Fig.[Tn|one can see the invariant yields of electrons from 
heavy-flavour decays in Au-Au and Pb-Pb collisions, for 



Table 3. The centrality classes and the corresponding average 
impact parameters at the kinematics of RHIC and LHC. 



Au-Au (^s = 


200 GeV) 


Pb-Pb (^ = 


5.5 TeV) 


C1-C2 


b (fm) 


C1-C2 


b (fm) 


0-10% 


3.27 


0-10% 


3.45 


10-20% 


5.78 


10-20% 


6.11 


20-40% 


8.12 


20-40% 


8.58 


40-60% 


10.51 


40-60% 


11.11 


60-92% 


12.80 


60-90% 


13.45 


0-92% 


8.44 


0-90% 


8.77 



a specific choice of the hydrodynamical scenario (viscous 
hydrodynamics with tq — 1 fm) and of the QCD scale fac- 
tor (/i — i'nT/2). Binary-scaled spectra are also shown 
as dashed lines. We are aware that cross sections are not 
the best observable for studying the effect of the nuclear 
medium on the propagation of the heavy quarks; for this 
purpose the nuclear modification factor is much more suit- 
able. It is however worth noting how the formalism is able 
to give a fair description of the PHENIX data over many 
orders of magnitude. Of course, this is also a consequence 
of the good performance of the pQCD outcome for the 
heavy quark production in NN collisions, as already dis- 
cussed in Sect 12.11 



3.2 Nuclear modification factor 

Let us consider now the nuclear modification factor, 

dNAA/dpT 



Raa{pt) 



Ncoi\dNpp/dpj 



(23) 



In order to explore the effect and the relative importance 
of those ingredients of the model that are not fully un- 
der control, that is the QCD scale factor and the hydro- 
dynamical scenario, it is instructive to start by studying 
the nuclear modification factor at quark level (R'aa) first, 
which is the one obtained from Eq. ([^^ by using the spec- 
tra before hadronization (keeping the same pseudorapidity 
window as for the electrons). Thus, in Fig. [Tl]we display 
the charm and bottom nuclear modification factors, for 
the case of viscous hydrodynamics (tq ~ 1 fm) and min- 
imum bias collisions, at both RHIC and LHC energies. 
Three choices for the scale factor n are considered: the 
strong dependence of Raa upon this parameter is quite 
apparent. We have not tried to fix the value of fi through 
a best fit to the RHIC data; however, in the following 
we employ the value fi — Z-kT /2, the one, among those 
displayed in Fig. [TTl that turns out to be in better agree- 
ment at large momenta with the PHENIX data on non- 
photonic electrons (see below). Another important effect 
on Raa can be appreciated by comparing the results at 
RHIC (panel (a)) and LHC (panel (b)) energies. Indeed, as 
already mentioned in Sect. 12. 1[ the elementary QQ pro- 
duction cross sections per binary collision entering into 
the numerator and denominator of Eq. (j23p . are not the 
same, because of nuclear effects. Raa for charm and bot- 
tom is then proportional to the ratio of the corresponding 
production cross sections. From Table [T] it clearly appears 
that for charm the relative play of shadowing and anti- 
shadowing is always giving rise to a quenching factor; for 
bottom, instead, one gets an enhancement at RHIC and 
a quenching at LHC (at ^/s = 5.5 TeV). 

In Fig. [T2] we display the charm and bottom nuclear 
modification factors for a fixed value ^ = 37rr/2 of the 
scale factor and several choices for the hydrodynamical 
scenarios (again for minimum bias collisions at RHIC and 
LHC). One can see that the uncertainty associated to the 
treatment of the hydrodynamical evolution of the QGP 
is quite modest compared to the one due to the heavy- 
quark interactions fFig. ITlT) . Specifically, the treatment of 
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Fig. 10. (a) Invariant yields of electrons [(e+ +e^)/2], from heavy-flavour decay in Au-Au collisions at ^/s = 200 GeV, at mid- 
rapidity as a function of transverse momentum (scaled for clarity). The experimental points represent PHENIX data [161117) : both 
statistical (error bars) and systematic (gray boxes) errors are displayed. The dashed curves represent the contributions coming 
from electrons originating from charm and bottom quarks generated by POWHEG with inclusion of transverse broadening and 
rescaled by A^coii. The solid curves are our results after propagation in the fireball (viscous hydrodynamics with tq = 1 fm, 
H = 3ttT/2). (b) Invariant yields of electrons [(e"*" -I- e~)/2], from heavy-flavour decay in Pb-Pb collisions at ^/s = 5.5 TeV, at 
mid-rapidity as a function of transverse momentum (scaled for clarity). Dashed and solid lines as in panel (a). 



the expanding fluid as either viscous or ideal has practi- 
cally no influence on the heavy-quark propagation. On the 
other hand, some sensitivity to the choice of the starting 
time of the hydrodynamical evolution is seen to exist. In 
fact, earlier times expose the heavy quarks to higher tem- 
peratures and give rise to a more pronounced quenching 
of Raa- At the LHC the ideal scenario with tq = 0.1 fm 
appears to be somehow in contradiction with the trend 
displayed by the other cases, showing smaller medium ef- 
fects for a shorter equilibration time. Note that it might 
just signal some shortcomings of the hydrodynamical code 
in a regime where it has received little testing. Up to semi- 
peripheral collisions the effect of an earlier equilibration 
time is anyway moderate and smaller not only than other 
theoretical uncertainties, but also than the experimental 
ones. However, as we shall see below, for peripheral col- 
lisions the effect due to the choice of tq turns out to be 
more important. 

In Fig. we display the charm and bottom nuclear 
modification factors of heavy quarks, open heavy-flavour 
hadrons and electrons for viscous hydrodynamics (tq = 
1 fm), scale factor /i = 3ttT/2 and minimum bias collisions 



at RHIC and LHC. The effect of hadronization and decay 
can be summed up as a general "softening" of the nuclear 
modification factor, especially after the decay of B mesons. 

In Fig. [T3] the results of our calculations of the nu- 
clear modification factor of open heavy-flavour electrons 
at RHIC (for the case of viscous hydrodynamics with tq = 
1 fm and /i = 37rT/2) are compared to the data from the 
PHENIX experiment [TOIITT] in all the centrality classes. 
Here we display the Raa's due to the electrons originat- 
ing either from c-quarks or from 6-quarks and their com- 
bination, showing also the statistical errors of the calcula- 
tions. For pt£3 GeV/c the data are generally underesti- 
mated: this is the region where, as already mentioned, the 
hadronization mechanism of coalescence has been shown 
to be important and to give rise to an enhancement of Raa 
[23] . At larger momenta the calculations turn out to be on 
the whole in agreement with the data from PHENIX, with 
an evident contribution from the bottom quarks. 

The above considerations are valid in all the centrality 
classes, except perhaps the most peripheral one (60-92%): 
here the effect of the medium is expected to be pretty 
weak and, indeed, our calculation does not show, at low 
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Fig. 11. (a) The heavy-quark nuclear modification factor at RHIC for viscous hydrodynamics (ro = 1 fm), h = 8.44 fm and 
different choices of the QCD scale factor; both c (light lines) and 6 (heavy lines) quarks are shown, (b) As in panel (a), but for 
LHC at 6 = 8.77 fm. 
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Fig. 12. (a) The heavy-quark nuclear modification factor at RHIC for /i = 37rr/2, b = 8.44 fm and different hydrodynamical 
scenarios; both c (light lines) and b (heavy lines) quarks are shown, (b) As in panel (a), but for LHC at 6 = 8.77 fm. 
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Fig. 13. (a) The nuclear modification factor of heavy quarks, open heavy-fiavour hadrons and electrons at RHIC for /i = 37rT/2, 
viscous hydrodynamics (ro = 1 fm) and b — 8.44 fm; both the charm (light lines) and bottom (heavy lines) sectors are shown, 
(b) As in panel (a), but for LHC at & = 8.77 fm. 
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Fig. 14. The nuclear modification factor of open heavy-flavour electrons at RHIC for ^ = 3nT/2 and viscous hydrodynamics 
{to — 1 fm) in various centrality classes. The circles are data from the PHENIX experiment [161117] . including both statistical 
(error bars) and systematic (grey boxes) errors; the other points (with the statistical errors shown) are the outcome of our 
calculations for electrons originating from charm and bottom quarks and their combination. 
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Fig. 15. (a) The nuclear modification factor of open heavy-flavour electrons from peripheral collisions at RHIC for fi = 3nT/2 
and different starting times for the hydrodynamical evolution. Data are from the PHENIX experiment [161117) . The theoretical 
curves show (without statistical errors) the total yield of electrons from charm and bottom quarks, (b) As in panel (a), but for 
minimum bias collisions. 



momenta, any underestimation of the data, but rather it 
overshoots them at intermediate momenta. In Fig. ITSl a) 
we compare the PHENIX data for the most peripheral 
centrahty class to the outcome of calculations based on 
different initial times. A reduction of tq from 1 fm to 
0.1 fm in all the other centrality classes gives rise to a 
rather mild effect, well within the theoretical and exper- 
imental uncertainties, as one can see, e. g., in Fig. I15f b) 
for minimum bias collisions. For very peripheral collisions, 
on the other hand, the highest temperatures associated to 
earlier equilibration times provides a stronger signature. 
Note, however, that for this centrality class — where one 
is moving from a region where the core of the nuclear den- 
sity distributions is probed to a region where only the tails 
of the density distributions are involved — the approxi- 
mation based on the use of an average impact parameter 
might be less reliable. 

In Fig. [16] we display the nuclear modification factor 
of open heavy-flavour electrons at LHC@5.5 TeV for QCD 
scale factor fi — 3ttT/2 and two choices of the equilibra- 
tion time. The features of Raa appear to be quite similar 
to the case of RHIC, except for a larger quenching at small 
momenta — due in our calculation to the smaller pp over 
AA ratio of the cc production cross sections — and for a 
slightly smaller value at large momenta. The choice of the 
equilibration time, as already discussed, has a mild effect, 
except perhaps on the most peripheral collisions. 

In Fig.[T7]the nuclear modification factor of open heavy- 
flavour electrons, obtained by integrating the electron yields 
above a given transverse momentum, is reported as a func- 
tion of the number of participants. When the integration 
range covers most of the spectrum, one should get the ratio 
of the total electron yields in AA and pp collisions, that is. 



in practice, the ratio of the respective heavy-quark produc- 
tion cross sections. When only large transverse momenta 
are selected, one observes the quenching of Raa due to the 
softening of the heavy-quark spectra by the medium. The 
calculations are again in fair agreement with the PHENIX 
data, except, as already noted, for the most peripheral col- 
lisions. These, however, represent a real puzzle, since they 
seem to imply that the strongest medium effects manifest 
themselves with the highest and the lowest numbers of 
participants. 

Finally, in this study we have also tried to understand 
the amount of uncertainty introduced by an approxima- 
tion that is, to our knowledge, always employed in this 
kind of calculations, i. e. the use of an average impact pa- 
rameter to represent a given centrality class. In Fig.lTSlthe 
minimum bias nuclear modification factor of open heavy- 
flavour electrons at RHIC and LHC is displayed, compar- 
ing the results obtained by using a single average impact 
parameter to those obtained by combining the weighted 
spectra in all the centrality (sub)classes (0-10%, 10-20%, 
...). The effect is not dramatic compared to other uncer- 
tainties, especially at LHC, but for an accurate estimate 
of Raa it will have to be accounted for. 



3.3 Elliptic flow coefficient 

The anisotropy parameter associated to the elliptic flow 
of the heavy quark is defined as 

fd0^/^cos20 
= , (24) 

'^'^ d(j)dpT 
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Fig. 16. The nuclear modiflcation factor of open heavy-flavour electrons at LHC@5.5 TeV for /x = 2>itT/2 and viscous hydro- 
dynamics (to = 1 fm and 0.1 fm) in various centrality classes. Electrons originating from charm and bottom quarks and their 
combination are shown, including the statistical errors. 
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Fig. 17. The nuclear modiflcation factor of open heavy-flavour electrons at RHIC (a) and LHC@5.5 TeV (b) for fi — 3nT/2 
and viscous hydrodynamics (tq = 1 fm), obtained by integrating the electron yields over the indicated momentum ranges, as 
a function of the number of participants. The data points from the PHENIX experiment [T7], including both statistical (error 
bars) and systematic (grey boxes) errors. 
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Fig. 18. The nuclear modiflcation factor of open heavy-flavour electrons at RHIC (a) and LHC@5.5 TeV (b) for /i — 3nT/2 
and viscous hydrodynamics (tq = 1 fm). The circles refer to the calculation at a flxed average impact parameter, the squares 
to the weighted sum of the results in all the centrality classes. 



where (PNAA/d4>dpT is the doubly differential yield of 
open heavy-flavour electrons, of D and B mesons or of 
c and b quarks, depending on which observable one is con- 
sidering. 

In Fig. [19] we display the charm and bottom ellip- 
tic flow coefficients of heavy quarks, open heavy-flavour 
hadrons and electrons for viscous hydrodynamics (tq — 
1 fm), scale factor fj, = 37rT/2 and minimum bias colli- 
sions at RHIC and LHC. In a momentum range up to a 
few GeV/c the heavy quark V2 is represented by a grow- 
ing function ofpr- Since, as already noted, hadronization 
and decay to electrons produce a softening of the respec- 
tive spectra, the net result is a stronger anisotropy of B /D 
mesons with respect to c/6 quarks and an even stronger 
one for the open heavy-flavour electrons. At the LHC en- 
ergy one gets larger elliptic flow coefficients, especially the 
bottom contribution. 

In Fig. [20] we display the total V2 of open heavy-flavour 
electrons for minimum bias events at RHIC and LHC. The 
calculation at RHIC energy is compared to the data from 



PHENIX, which, however, are available only up to pt ~ 
4 GcV/c. Our results for tq = 1 fm underestimate the data 
for Pt < 3 GeV/c, as they also do in the case of the nuclear 
modiflcation factor. The already mentioned coalescence 
mechanism for the hadronization process is known to be 
important in this momentum region and to provide an 
enhancement of both Raa and V2 [M]. Note, however, that 
an earlier equilibration time gives rise to an anisotropy 
parameter closer to the data. At the LHC energy, one 
observes a stronger anisotropy for the exclusive c and b 
contributions, while the effect due to the choice of tq turns 
out to be milder. 



4 Summary and conclusions 

In this work we intended to provide a consistent study of 
the stochastic dynamics of heavy quarks (c and b) in the 
hot environment formed in high-energy heavy-ion collision 
experiments. The background medium has been taken as 
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Fig. 19. (a) The minimum bias anisotropy parameter V2 of heavy quarks, open heavy-flavour hadrons and electrons at RHIC 
for n = 37rT/2, viscous hydrodynamics (ro = 1 fm) and b = 8.44 fm; both the charm (light lines) and bottom (heavy lines) 
sectors are shown, (b) As in panel (a), but for LHC at fe = 8.77 fm. 
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Fig. 20. (a) The minimum bias anisotropy parameter V2 for open heavy-flavour electrons at RHIC for fi — 3ttT/2 and viscous 
hydrodynamics (tq = 1 fm and 0.1 fm). The circles are data from the PHENIX experiment [161117) . including both statistical 
(error bars) and systematic (grey boxes) errors; the other points (with the statistical errors shown) are the outcome of our 
calculations for electrons originating from charm and bottom quarks and their combination, (b) As in panel (a), but for LHC 
and two choices of ro . 



a fluid (the QGP locally in thermal equilibrium), whose 
evolution is governed by relativistic hydrodynamics, while 
the interaction of the heavy quarks with the plasma has 
been treated within a weakly-coupled framework. 

To accomplish this programme one has to evaluate the 
relevant transport coefficients (accounting, in our case, for 
the squared momentum acquired per unit time) charac- 
terizing the dynamics of the heavy quarks coupled to the 
medium. We have done this by summing the contributions 
of soft and hard collisions in the plasma. The first ones 
have been described in the HTL approximation, while the 



latter have been treated within a kinetic pQCD calcula- 
tion. The propagation of the heavy quarks in the expand- 
ing plasma has been described by a relativistic generaliza- 
tion of the Langevin equation. 

When the hydrodynamical evolution of the QGP has 
reached the region of energy density where the quark to 
hadron transition may occur, the heavy quarks are let 
hadronize and eventually decay into electrons: in fact, 
their spectra at RHIC are so far the only source of infor- 
mation on the heavy-flavour dynamics in the QGP (while 
waiting for the forthcoming open-charm analyses at RHIC 
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and LHC). Results have been presented for the invariant 
single-electron spectra, for the nuclear modification factor 
Raa and for the elliptic flow coefficient V2- 

The main goal of our work has been first of all to 
provide of a fully consistent weakly-coupled calculation 
(accounting for medium effects at the lowest non-trivial 
order), to be viewed as a benchmark for more advanced 
studies or less conventional scenarios. Nevertheless, the 
overall degree of agreement with the current experimen- 
tal data looks quite satisfactory, suggesting that the ac- 
tual values of the heavy-quark transport coefficients might 
be not too different from what predicted by a weak cou- 
pling approach. In particular, for large enough values of 
Pt (say pt^3 GeV/c), it has been possible to reproduce 
the experimental quenching of the single-electron spectra 
for any centrality class (except perhaps the most periph- 
eral events, which will be discussed in more detail in the 
following). On the other hand, at low pT (<;3 GeV/c) our 
results for Raa are consistently below the data, predict- 
ing too much quenching. This outcome might be a conse- 
quence of the implemented hadronization scheme, where 
we have accounted only for the fragmentation mechanism 
(which entails a degradation of the heavy-quark momen- 
tum). Hadronization through coalescence with light par- 
tons of the medium could help in getting closer to the 
experimental data. Hadronization represents therefore an 
important source of systematic uncertainty for moderate 
values of pr ■ 

Concerning the sensitivity to the properties of the back- 
ground medium, in general we have found a rather mild 
dependence on the hydrodynamical scenario (ideal or vis- 
cous, with different choices of tq) used to describe the 
plasma, except for the most peripheral events, which de- 
serve some separate comments. Quite surprisingly, though 
with large error bars, experimental data on heavy-flavour 
electrons from RHIC display a sizable quenching of the pr- 
spectra even in the 60 — 92% centrality class. As shown in 
Sect. 13. 2[ employing a value of tq ^ 1 fm, as often done 
in hydrodynamical simulations at RHIC energy, would 
lead to Raa ~ 1, thus essentially showing the absence 
of medium effects. On the other hand, assuming a much 
faster thermalization, by setting tq = 0.1 fm, it is possi- 
ble to get much closer to the data. Notice that with such 
a small value for the equilibration time it is also possi- 
ble to reach a slightly more satisfactory agreement with 
the elliptic flow V2 measured for heavy-fiavour electrons at 
RHIC. At the present stage we cannot draw any definite 
conclusion on this point and we leave it as an open issue 
for future investigation. 

Besides analyzing the RHIC data, we have also pro- 
vided predictions for LHC (at the highest energy that 
should be reached, \/sj^j^ — 5.5 TeV). In the absence of 
experimental data to constrain the properties of the back- 
ground medium, we have explored some of the possible 
hydrodynamical scenarios proposed in the literature. As 
a general outcome, charm is found to display a stronger 
quenching and a much larger elliptic flow with respect to 
RHIC. Also the spectra of bottom are more suppressed 
and characterized by a modest elliptic flow. However, in 



the combined electron spectra (ec+f,), due to the larger 
relative weight of the b contribution, differences between 
RHIC and LHC are less evident. The possibility by the 
ongoing experiments of identifying the separate contribu- 
tion of charm, by looking at its hadronic decay channels, 
will certainly increase the amount of information on the 
heavy-fiavour dynamics provided by the data. Heavy-ion 
collisions at LHC are presently being performed at the 
NN center-of-mass energy \/s^jv — 2. 76 TeV: work is in 
progress in order to analyze, within the framework pre- 
sented in this article, also this kinematical setup. 

We thank R. Averbeck, G. Ridolfl, A. Toia and J. Bielcik for 
helpful discussions and P.B. Gossiaux for providing us his data 
on the transport coefficients. 



A Transverse momentum broadening of 
heavy quarks in nucleus-nucleus collisions 

As explained in the text, our procedure to account for the 
transverse momentum broadening in a collinear collision 
consists in adding, to the heavy quarks generated in each 
event, a transverse component randomly chosen according 
to a Gaussian distribution (41j 

5(M = ^^exp(-4/(4)). (25) 

As already mentioned, for pp collision the typical value of 
(fc2)NN is^l GeVVc^ 

Let us now consider a collision between the nuclei A 
and B at impact parameter b, in which a QQ pair is pro- 
duced in a hard event occurring at position s in the trans- 
verse plane (measured with respect to the center of nu- 
cleus A) and at the longitudinal coordinates (measured in 
the rest frames of the respective nuclei) za and zb- The 
heavy-quark pair will be produced by partons (mainly glu- 
ons) that have already crossed a certain thickness of nu- 
clear matter, thus having acquired some transverse mo- 
mentum. In Refs. [44[|45| the average squared-momentum 
of the produced QQ heavy meson has been found to be 
given, neglecting absorption effects, by the expression 

(p|)2f(b,s,z^,ZB) = (p|)g^ (26) 
+agN[lA{s,ZA) + Ib{s - b, Zb)], 

where {p't)nn ^^'^ value for the same quantity in NN 
collisions and 

dz pa{s,z)/po 

-oo 

/• + 00 

Ib{s -b,Zb) = j dzpB{s-b,z)/po (27) 

J Zb 

represent the lengths traveled by a gluon of B in nucleus 
A and vice versa. In Eq. (|27)) p{s, z) is the nuclear density 
function and po is the central nuclear density. The pa- 
rameter agM represents the average squared (transverse) 
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momentum per unit length acquired by a gluon in nuclear 
matter. At the energy of SPS one finds good agreement 
with the J/i: data for agN « 0.08 GeVVfm gS]. 

Let us now apply the above scheme, developed for char- 
monia, to the problem of inclusive single-particle spectra. 
The random transverse momentum kick of the pair re- 
flects the one acquired by the incoming gluons in crossing 
nuclear matter. It is shared by the two quarks, 



Table 4. The total nucleon-nucleon inelastic cross section and 
the parameter a^^ at different kinematics. 



^ (GeV) 


158 200 5500 


CTNN (mb) 


33 42 72 


a-^ (GeVVfm) 


0.081 0.10 0.17 


a^'^ (GeVVfm) 


0.221 0.27 0.47 



(28) 



and for a random choice of the azimuthal angle of each 
quark of the pair one has: 



1\Q 



COS 



= 2{kl 



(29) 

In general one does not know the longitudinal coordinate 
of the hard event, so that one takes an average over the 
longitudinal positions za and zb, with a weight given 
by the corresponding local nuclear density. The average 
squared momentum used in the distribution of Eq. (j25p 
to extract the transverse momentum kick for the heavy 
quarks is then: 



'"/-c^ '^^A Pa{s, za)Ia{s, za) 

I-^^^A PAis,ZA) 



a-gN 

2 

J^^ dzB Pb{s - b, zb)Ib{s ~ b, zb) 
f^^ dzB PB{s-b,ZB) 



(30) 



Note that we keep track of the transverse position of the 
hard event, since the transverse momentum kick is im- 
parted to heavy quarks distributed in the transverse plane 
according to binary-collision scaling. 

The parameter agM has been fixed in Refs. [^H5] at 
an energy lower that the ones we are interested in (-^s = 
158 GeV gS]). In order to extrapolate its value to the 
RHIC and LHC regimes we have considered the transverse 
squared momentum of Eq. ([50)) averaged over the QQ pair 
transverse position s and over the impact parameter b. 
This same quantity has been modeled in Ref. .62u41| in 
terms of the squared-monientum transfer per collision. 



A'^ip) = 0.225 



ln^(A^/GeV) 
l + ln(^/GcV) 



GeV' 



(31) 



(/i = 2mQ being the interaction scale) and of the inelas- 
tic nucleon-nucleon cross section unn- A straightforward 
comparison of our expression with the one of Ref. (41) 
yields, indeed, 

o-gN — ^^ctnnPo, (32) 

(Tnn being the total inelastic NN cross section and po the 
central nuclear density, and leads to the values displayed 
in Table H 



B Momentum difFusion coefficients: 
calculation of the hard contribution 

In this Appendix we provide some details on the calcula- 
tion of the hard-collision contribution to the momentum- 
diffusion coefficients. The squared amplitude for scattering 
on quarks (Fig. |4]) reads: 



(33) 



{NP - uf + {s- M^f + 2MH 

where, as mentioned in the text, the factor 2A^/ arises 
from the equal contribution of all light quark and anti- 
quark flavours, while the further factor 2A^c with respect 
to Ref. ^29 reflects our summing, rather than averaging, 
over the helicities and colours of the incoming light quark. 
For scattering on gluons (Fig. 14]) one gets 

\Mg{s,t)\^ ^2{Nl~l)al^^ 
is - ]VP){AP - u) 



32- 



t2 



64 (s - hP){KP - u) + 2M2(s + hp) 

~^ (S - M2)2 

64 (s - M^){]VP - u) + 2]VP{AP + u) 
16 



(AP 

]VP{A]VP - t) 



m)2 



9 (s-M2)(M2-u) 



+16 
-16 



(s - AP){AP -u) + KP{s - u) 

t{s-AP) 
(s - Af2)(M2 -u)- M2(s - u) 



t{A'P - u) 



(34) 



where the factor 2(iV2 — 1) with respect to Ref. [53] stems 
from summing over the polarization and the colours of the 
incoming gluon. 

In order to perform the phase space integration in 
Eqs. (IT5|) and p^ . it is convenient to express the lon- 
gitudinal and transverse squared momenta exchanged in 
the collision in terms of the energy uj lost by the heavy 
quark and of the Mandelstam variable t. One has 



{2Eu:+\t\f 
4p2 



and 



+ \t\ 



{2Euj+\t\)' 



4p2 



(35) 
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SO that one can follow the procedure of Ref. [12] for the 
evaluation of integrals like (J^, = J dfc/(27r)^) 



nn 



1 fnB/F{k) f l±nB/F{k') f e{\t\ - \t\*) 



2E 



2k 



2k' 



p' 



2E' 



X (27r)4(5('^) {P + K-P' - K')f{s, t, w). 



(36) 



One gets 

m = 

where 
being 



nB/F{k) /"I*! 



Wn'^pE Jk 2k 7| 



d\t\ 



did 



x[l±nB/F{k + uj)]f{s,t,u}), 

g{uj) — —a?uP' + 5w + c, 



(37) 
(38) 



a=^^, h='^[E{s-M^)~k{s + M^)\,m 



P 



\t\ 



c = U [4p2;(,2 _ _ j^j2 _ 2Ek) - \t\{{E + kf - s)] 



and 



^inax/min 



2^^ ■ 



From + Aa^c > one gets then: 



1^1 



2\2 



(40) 



(41) 



Hard collisions provide a non- vanishing contribution only 
as long as > This allows to set the relevant 

range of integration over fe, which is given by: 



k > 



1 + Wi 



UP 



(42) 



1 1 + ^ 

'v AEkv I V 1*1* 



AE{l + v) 
— 1 < cos Okp < min 
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